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Q\ ' The excess of high-Q 2 events found by HI and ZEUS at HERA in e + p deep-inelastic 

scattering above the Standard Model prediction motivates us to calculate the NLO 

O-i! 

QCD corrections to the HERA scalar leptoquark (or squark) production cross see- 
On. 

tions. We find that the corrections are significant, of order 50% in the mass range 



of interest. We also calculate the leptoquark average transverse momentum squared 
and find it to be rather small. Various leptoquark production cross sections at the 



Tevatron pp collider are also considered. We investigate in detail the leptoquark 
interpretation of the HERA data. First we assume a minimal leptoquark model 
with a single diagonal Yukawa coupling to first family lepton and quark mass eigen- 
states only. In this case constraints from atomic parity violating experiments allow 
only isodoublet scalar leptoquark production at HERA. This interpretation can be 
confirmed or ruled out in the near future by high luminosity data at the Tevatron. 
The Tevatron data already appear to rule out the vector leptoquark interpretation 
of the HERA data. We also consider a more general model which allows for all pos- 
sible left-handed, right-handed, flavour and lepton number changing couplings. The 
allowed values of the Yukawa couplings of this general model offer several different 
interpretations of the data which are radically different from the minimal model 
solutions. However these somewhat ad hoc tuned solutions can easily be tested by 
future HERA experiments. 



1 E-mails: kunszt@itp.phys.ethz.ch; W.J.Stirling@durham.ac.uk. 



1 Introduction 



Recently the HERA collaborations HI |l] and ZEUS f|| have presented new results on 
high-Q 2 deep inelastic e + p scattering based on 14 pb~ x and 20 pb _1 luminosity data, 
respectively. The data indicate an excess of events with respect to the Standard Model 
expectation predicted by the conventional next-to-leading order QCD-improved parton 
model in a new kinematical range unexplored by previous experiments. At high Q 2 values 
Q 2 > 15000 GeV 2 HI finds 7 excess events while ZEUS finds 5 events. An interesting 
feature of the HI events is that they are clustered in a single bin of the invariant mass 
of the jet-lepton system (M e j = yfxs) around 200 GeV with an estimated experimental 
resolution of 6 GeV The ZEUS events do not exhibit such a clustering, but this is 
not in conflict with the HI data since ZEUS have a larger M e j resolution. HI have also 
published results on charged current events |L|. Although they find more events than 
expected, the excess is not statistically significant. If the signal is eventually confirmed as 
the production of a new heavy resonance, one will presumably be able to measure its spin 
and its branching ratios both to charged current and to flavour and/or lepton number 
violating channels. 

The possibility of s-channel resonance production at HERA has been suggested in two 
theoretical schemes. First, Buchmuller and Wyler || argued that the leptoquarks (LQ) 
necessarily appearing in grand unified theories [HI [5| can also be accommodated in theories 
with conserved lepton and baryon number. They pointed out that these leptoquarks 
can be light (i.e. 771lq <C Mqut) since low energy experiments impose much weaker 
bounds on lepton and baryon number conserving couplings. Second, the squarks (q) 
of supersymmetric extensions of the Standard Model can have direct (lepton number 
violating) couplings to quarks and leptons if R-parity is not conserved. Therefore they 
can also be coupled to lepton-quark initial states and thus appear as s-channel resonances 
0. In the present study we shall restrict our discussion to the leptoquark interpretation 
only, although several of our results will remain valid also for squark production. 

In section 2 we consider leptoquark production cross sections at HERA, including next- 
to-leading order QCD corrections. From the number of excess events we extract the value 
of the leptoquark couplings (A) and discuss the constraints on these coming from atomic 
parity violating experiments. We also comment on non-minimal models which allow for 
all possible Yukawa couplings, following the analysis of ref. 0. 

In section 3 we consider single and pair production of scalar and vector leptoquarks 
in pp collisions at the Tevatron. We study A dependent contributions, the effects of QCD 
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corrections, signals and backgrounds, and experimental limits on the leptoquark masses 
obtained by DO and CDF. 

Section 4 describes the calculation of the next-to-leading order QCD corrections to 
colour-triplet scalar resonance production in ep collisions. The numerical value of the 
'i^-factor' and the average transverse momentum squared of the produced resonance are 
also studied. Full details of the calculations are given in the Appendix. Our conclusions 
are contained in section 5. 

2 Leptoquark production at HERA and limits on 
couplings 

We first consider the minimal scheme in which the leptoquarks have separate baryon 
and lepton number conserving chiral couplings for each family of mass eigenstates. This 
allows us to evade severe low energy limits on lepton number violating couplings J7|, ||, 
[J. Leptoquarks which are coupled to the first family only are called first generation 
leptoquarks. A leptoquark is said to couple chirally if it couples either to left-handed (L) 
or right-handed (R) leptons but not both. 

The allowed leptoquark representations have been classified in ref. || . There are seven 
B and L conserving and SU(3) xSU(2) xU(l) renormalizable couplings both for scalar 
leptoquarks 

£s = giL^L^hSi + giRU c R ir 2 e R S 1 + giR^R^enSt + g 2R u R l L S 2 (1) 
+ g2Lq L e L S 2 + g2L~d R l L S 2 + g^LQl^hS^ + h.c. , 
and vector leptoquarks 

£v = hi L q L ^l L V^ + h 1R u R 'j^e R Vl 1 + h 1R u R 'j^e R Vl 1 + h 2R d R i^ l J, L V^ (2) 
+ hiL^hl^By^ + ^L^Rlixhy% + h iL q L -i^fl L V^ + h.c. . 

We can distinguish the leptoquarks according to their weak isospin properties: singlets 
(e.g. Si), doublets (S 2 ) and triplets (S3). The doublet scalar leptoquarks and the sin- 
glet and triplet vector leptoquarks have zero fermion number (L + 3B) while the singlet 
and triplet scalar leptoquark and doublet vector leptoquark have fermion number two. 
Couplings are designated left or right handed according to the handedness of the lepton. 
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We emphasize a difference in the theoretical status of scalar and vector leptoquarks. 
Within the Standard Model, vector leptoquarks cannot be considered as gauge bosons and 
therefore their couplings to the gauge bosons of the Standard Model are not renormaliz- 
able in general. Scalar leptoquarks, however, can be considered as possible renormalizable 
matter fields of the Standard Model. In particular, the squarks of the supersymmetric 
standard model with R-parity violation provide a theoretically appealing realization of 
the isospin doublet and singlet leptoquarks (S 2 and Si). In contrast, scalar leptoquarks of 
exotic electric charge (Si), and S3 and vector leptoquarks have no counterparts in super- 
symmetric models. In section 3 it will be shown that the Tevatron limits on the production 
of vector leptoquarks are only very marginally consistent with the vector leptoquark in- 
terpretation of the HERA events. In this section, therefore, we shall consider limits and 
production rates only for scalar leptoquarks. 

The Yukawa couplings in eq. (|l|) are diagonal in the family number. In order to avoid 
the severe low energy limits from the leptonic decays of charged pions it is convenient to 
require that the leptoquarks have only one non-zero chiral coupling. Thus we can have 
left handed or right handed scalar leptoquarks Sir, Su, and S 2 l,S 2 _r in flU). In Table 1 
we list all possible scalar leptoquark states which can couple to a e + q initial state. 

The leading order production cross section at HERA has the simple form (see Appendix 
A.l) 

a = ^- q (ml Q / s ) ■ BR , (3) 

where A 2 denotes generically the coupling constants squared {g\ R ) of the produced scalar 
leptoquark as listed in Table 1 (or two times the coupling constants squared [h\ R ) of the 
vector leptoquarks). The cross sections of scalar leptoquark production calculated to next- 
to-leading order (0(cns)) accuracy are plotted in Fig. 1 as a function of the leptoquark 
mass ttt-lq for the e + u, e + d, e + s, e + u and e + d initial states^ assuming A 2 = 1. For other 
values of A the cross section obviously scales as A 2 . The cross section values vary strongly 
according to the initial quark flavour. For ttilq = 200 GeV we find the ratios (see section 4, 
Table 3) 

a u : a d : o~ d :a n :a s = 1440 : 352 : 16.7 : 7.32 : 8.21 (4) 

where the values are given in picobarns. Assuming 80% experimental efficiency, we can 
extract from the HI and ZEUS results a rough estimate for the HERA leptoquark pro- 
duction cross section of 0.7 pb. Comparing this value with the predictions of Fig. 1 for 

2 At 0(as) the e + g scattering process also contributes, see Appendix, and is included in the cross 
sections in Fig. 1. 
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Table 1: Scenarios for scalar leptoquark production in e + p collisions, with values and 
limits for the various couplings obtained from HERA and atomic parity violation (APV). 
Limiting values for the branching ratios corresponding to definite chiral couplings (i.e. L 
or R) are given. 



m LQ = 200 GeV we can calculate numerical values for the various coupling constants listed 
in Table 1, assuming branching ratios which correspond to definite chiral couplings. 

In the the minimal scheme the strict bounds on electron and muon number conservation 
are automatically satisfied. Allowing only one chiral coupling, the severe constraints from 
the leptonic decays of charged pions and kaons are also avoided. But in the case of 
first generation leptoquarks the diagonal flavour conserving couplings (see Table 1) are 
also severely limited by the precise measurements of the atomic parity violating neutral 
current weak charge of the caesium atom 

Q w = -376C lu -A22C ld , (5) 

where C\ u and C\d denote the lepton axial vector - quark vector interference terms of the 
effective low energy four fermion interaction 

£apv = 4 E Cuh^lq^Qi, (6) 

with 

C lu = -- + -sin 2 6 w , C ld = - - - sin 2 Q w . (7) 
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The scalar leptoquark contribution to C\ q is 



V2g 



2 



ACl « SG F ml Q 



The measured value is 10 



Q w = -71.04 ± 1.58 ± [0.88] th . (9) 
to be compared with the Standard Model prediction ]10 



Q w = -72.88 ± 0.05 ± 0.03. (10) 
The difference between theory and experiment is therefore smaller than 

|AQ W | = 3.6, (11) 

which gives 

° 2 <°H2^k) 2 < 12 > 

which in turn leads to the values given in the fourth column of Table 1. From this we 
conclude that the HERA and the APV data are consistent only with the production of 
isodoublet (zero fermion number) leptoquarks. The parton distributions for antiquarks 
are too small at large x (~ 1/2) to allow sizeable production of leptoquarks with the 
weak coupling required by atomic parity violation data. If this interpretation is correct 
and only one leptoquark exists in the HERA kinematic range, then it follows that in 
e~p collisions the signal will be weaker by about two orders of magnitude! While this is 
an attractive explanation by itself, it may be excluded soon by new analyses from CDF 
and DO at the Tevatron. We shall see in the next section that scalar leptoquarks of 
mass mLQ ~ 200 GeV are very close to the published discovery limit for leptoquark pair 
production at the Tevatron. 

The Tevatron limit can be weakened by suppressing the branching ratio into charged 
lepton + jet, since the production is essentially independent of the leptoquark coupling. 
This motivates us to consider the more general parametrization of leptoquark production 
as studied by Davidson et al. J7| , with general mixed couplings between the three families. 
In the Lagrangian £,$ of (|TJ) the family diagonal leptoquark couplings g are replaced by the 
coupling matrix g % i where the summation labels run over the family indices i,j = 1,2,3. 
Davidson et al. summarize their results^ in terms of bounds on coupling constant 
3 See Table 15 of ref. @. 
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Table 2: Generation-mixing leptoquark coupling limits, from ref. J7|. The interactions 
denoted by * involve a top quark and the corresponding constraints have to be evaluated 
using the known top quark mass. 



combinations 



\gijgkm ( "'l-Q 1 



2" a VlOOGeV, 
In Table 2 we reproduce their results for g l i g mn at rn-^Q = 200 GeV. 

We see from Table 2 that there are many ways to avoid the bounds and still maintain 
consistency with the HERA signal. It is still premature to make a detailed study, therefore 
we mention only one possibility. Suppose that all the couplings are extremely small except 
(12) (12), i.e. when the leptoquark is coupled to a strange quark and a positron or to a 
charm quark and a positron. The smaller strange quark and charm quark content of 
the proton requires larger leptoquark couplings (see Fig. 1) but atomic parity violation 
constraints are avoided since there is no coupling to a second generation quark in the 
proton. In this scenario the signal in e~p scattering will be as large as for e + p. If we 
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insist on a smaller branching ratio (< 1) then we may allow either for both left and right 
handed couplings or for a significant (12) (23) coupling. Another possibility is to allow 
also for a (12) (31) coupling. In this way the Tevatron mass limit could be relaxed. In this 
scenario charged current decay modes can also easily be arranged. Finally we note that 
(i) the (12) (12) solution resembles the squark solution suggested in |IT], and (ii) although 
the above choice may appear somewhat ad hoc we should recall the huge hierarchy of the 
Yukawa couplings in the Standard Model without scalar leptoquarks. 



3 Leptoquark production at the Tevatron 



Colour-singlet leptoquarks have a standard gauge coupling (g s ) to gluons and can therefore 
be copiously pair produced in hadron-hadron collisions. The direct coupling to quarks (A) 
gives additional contributions to the production cross section, but these are numerically 
much less important for the range of A allowed by the HERA and APV data. In addition, 
vector leptoquarks can have additional anomalous couplings (/ty, Ay) to gluons, analogous 
to the anomalous electric quadrupole and magnetic dipole moments of W bosons. Al- 
though these lead to a violation of unitarity in the production cross section, they may 
be regarded as an effective low energy parametrization of a more complicated theoretical 
structure. 

Assuming that LQ — > I + q is the dominant leptoquark decay mode, the signature 
of pair production in hadron collisions - two energetic leptons and two energetic jets 
widely separated in phase space - is rather distinctive. Backgrounds from processes like 
W, Z+jets, 66+jets, etc. can be suppressed in principle by kinematic cuts. 

The pair production subprocess cross sections are 4 , with [3 2 = 1 
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for scalar pair production and |13|, [14|, 115, 16[ |17 



a 99->VV 



TCCXc 



96m2 Q 



3 (f - + ^ 



— (65 - 83/? 2 + 19/? 4 - /3 ( 



In 



4 The cross sections for scalar pair production have been derived in the literature in many different 
guises. See, for example, ref. and references therein. 



7 



a qq^VV 



Ana 2 (P_ 
24 



23 - 3p 2 + 



(14) 



l-(3 2 _ 

for vector pair production, with Xy = Ky = 0. Note that in general the scattering am- 
plitudes for qq — > SS, VV receive additional contributions (cr ~ ctsA 2 , A 4 ) from diagrams 
involving direct quark-leptoquark couplings (e.g. qq — > LQLQ by i-channel exchange of 
a lepton). However for A 2 = (9(10 _1 ) or smaller, these contributions are much smaller 
than the 0(a 2 s ) contributions listed above. The generalization of (0) including non-zero 
anomalous couplings is presented in ref. [ |L3] . 

Fig. |3| shows both cross sections at = 1.8 TeV as a function of the leptoquark 
mass. Note that for ttt-lq ~ 200 GeV the qq annihilation subprocesses are dominant. The 
parton distributions are the MRS(R2) set from ref. (with as(M%) = 0.12) and the 
renormalization and factorization scales are set equal to wilq. We see from Fig. ||] that 
the cross sections for tti^q ~ 200 GeV vector and scalar pair production are O(10) pb 
and O(0.1) pb respectively. The former are therefore ruled out by the current Tevatron 
data |12[. The only caveat to this is that it is possible to 'fine-tune' the anomalous 
couplings Ky, Xy to suppress the vector pair cross section by between one and two orders 
of magnitude, see for example ref. [13]. However this would appear to be completely 



unnatural. The cross sections in Fig. |3] do not include NLO QCD corrections. Estimates 
for these can be extracted from the calculation of ref. |]I9| for squark pair production in 
the infinite gluino mass limit. Numerically, the NLO corrections increase the cross section 
by a modest O(+10%) in the mass range of interest (see Fig. 19(b) of ref. [jl9|1 ) . 

A 200 GeV scalar leptoquark would therefore give rise to approximately 10 events in 
a sample corresponding to 100 pb -1 at the Tevatron. In fact a recent detailed analysis by 
the DO collaboration [|T^] of their combined Run I data gives a (95%cl) lower mass limit of 
175 GeV (147 GeV) assuming BR(LQ -> eg) = 1(0.5).Q Thus the m LQ ~ 200 GeV (scalar 
LQ) interpretation of the HERA excess events is allowed, by a small margin. Obviously 
the Tevatron limit is weakened if the branching ratio to eq is lowered (see section 2 above). 

Finally, we show also in Fig. [3] single scalar leptoquark production cross sections from 
the subprocess qg — > LQ + / p0|: 



cr, 



48s 



1 + 6x — 7x 2 + 4x(l + x)\nx 



(15) 



with x = m^n/s < 1. In Fig. ^|we have taken, for illustration, A 2 = 0.01 and distinguished 

the cases where the leptoquark couples to u, d and s quarks. Thus, for example, for a 

5 There is as yet no combined Run I published first generation leptoquark mass limit from the CDF 
collaboration. 
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scalar leptoquark which is produced at HERA via e + u and has a 100% branching ratio 
into the same state, the curve labeled 'SJ gives the number of gu — > (e + u) + e~ and 
gu — ► (e"u) + e + events. Other A 2 values and branching ratios can be obtained by a 
simple rescaling of the curves. Note also that for a scalar leptoquark which gives rise also 
to charged current events are HERA, for example via e + d — ► S — > P e tt, the additional 
'missing energy' processes gd — > (P e iz) + e~ and git — > (v e u) + u e are possible. The 
corresponding event rates are readily estimated from the curves in Fig. |3|. 

The most important conclusion from the single leptoquark production curves in Fig. |3| 
is that for m-^Q ~ 200 GeV, the single production cross section is less than the pair 
production cross section for A 2 < (9(0.07, 0.16, 0.8) for coupling to u, d, s quarks. 



4 Production cross section in NLO at HERA 

In order to get an accurate value of the leptoquark coupling constant from the HERA 
data one should take into account the next-to-leading order QCD corrections. This is also 
required for consistency if parton distributions obtained with a NLO fitting procedure are 
used. In this section we describe the calculation for scalar leptoquarks since, as we have 
seen, vector particles are disfavoured by searches at the Tevatron. Note that the transverse 
momentum (kx) distribution of the leptoquarks with respect to the beam direction is a 
by-product of the full NLO correction. 

The NLO-corrected cross section can be written as 
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(16) 



where x = ttv^q/s and Cp = 4/3, Tr = 1/2. The coefficient functions are calculated from 
the 0(cts) Feynman diagrams. Using MS factorization and renormalization we obtain (see 
Appendix for more details) 
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(18) 
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Here /Mj R and fifjy denote the factorization and (ultraviolet) renormalizationP] scales re- 
spectively. In our numerical calculations we will use 



= Ar = PUV = m LQ • ( 19 ) 

Note that the QCD correction is universal for the different types of scalar leptoquarks 
discussed in section 1. 

Figure [I] shows the NLO-corrected production cross sections for e + q — > LQ with 
q — u,d, u, d, s at y/s = 300 GeV. For ease of comparison we have set the overall leptoquark 
coupling to unity, i.e. A 2 = 1. The parton distributions are the MRS(R2) set from ref. |fT8 



with a 5 (M|) = 0.120. The cross section hierarchy mainly reflects the hierarchy of the 
quark distributions at large x. The ratio of the NLO to LO contributions (the 'K-factor' 
JC) is non-negligible. For the valence quark cross sections, K increases from about 1.25 
to 1.45 in the mass range of interest, 175 GeV < /tj-lq < 225 GeV, while for the sea 
quarks K increases from about 1.3 to 1.6 in the same mass range. Table ^ lists cross 
section values for leptoquark masses relevant to the HERA high Q 2 events. In the limit 
^lq — > y/s, i.e. x — > 1, the correction is dominated by a soft-gluon double logarithm, 
K ~ 1 + a s C F \n 2 (l -x)/n. 





a(e + u) 


a(e + d) 


a(e + u) 


a(e + d) 


<j(e + s) 




= 175 GeV 


3289 


1026 


47.0 


99.6 


49.7 






(2659) 


(801) 


(35.7) 


(73.8) 


(37.9) 




= 200 GeV 


1436 


352 


7.32 


16.7 
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(337) 


(61.3) 
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(0.479) 



Table 3: Leptoquark cross sections (in pb) at NLO in e + p collisions at yfs = 300 GeV, 
assuming unit coupling A 2 = 1. The numbers in brackets are the corresponding leading- 
order cross sections. 



The transverse momentum distribution da jdk\ and the average transverse momentum 

squared (k 2 -,) of the leptoquarks can also be calculated from the 0(as) diagrams for the 

6 The leptoquark coupling A 2 is rcnormalized by the 0(as) vertex corrections, see Appendix, hence 
A 2 = X 2 (fj,i /V ). 
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processes eq — > gLQ and eg — > gLQ (see Appendix A. 4). For the average we obtain 
(4) vrA 2 a s (li 2 ) f 1 dz 
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Asa 2ix 



x Z 
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(20) 



where a is defined in eq. (|i"6|) and 
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Figure ^ shows y (kj.) as a function of rn-^Q for the various e + q production mechanisms. 
Note that the average transverse momentum is in general very small, because the emission 
of an additional energetic parton is heavily suppressed. 



5 Conclusions 

We have calculated the NLO QCD corrections to the production cross-section of a coloured 
s-channel scalar resonance (leptoquark or squark) at HERA and reevaluated various lep- 
toquark production cross sections for Tevatron. The average transverse momentum of the 
leptoquark produced at HERA was also calculated from the NLO result. Our calculation 
was motivated by the leptoquark interpretation of the excess of events found by HI and 
ZEUS at HERA in e + p deep-inelastic scattering above the Standard Model prediction. 
The simplest and perhaps most attractive 'new physics' interpretation of the data is in 
terms of a first family leptoquark with a single Yukawa coupling. Constraints from data 
on atomic parity violation, double beta decay and various flavour changing processes are 
consistent with the HERA data only for isodoublet scalar leptoquark production. This 
interpretation, however, is seriously challenged by severe constraints from leptoquark pair 
production at the Tevatron pp collider which is largely independent from the Yukawa cou- 
pling of the leptoquark. The Tevatron limits are even more severe for vector leptoquarks. 
Because the QCD corrections to the HERA cross sections are large and positive, the ex- 
tracted couplings are somewhat smaller than those obtained from a leading order analysis 
only. 

If one also considers leptoquarks which couple to several generations, then in the 
enlarged space of couplings one can find several solutions compatible with the low energy 
constraints and the HERA excess events. In particular, one can have solutions with 
reduced branching ratios and/or charged current events, as well as lepton number violating 
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processes. Furthermore one can no longer rule out the isosinglet and isotriplet scalar 
leptoquark interpretation in such models. The solutions with reduced branching ratios 
also relax the Tevatron limits. However, they do require a rather ad hoc hierarchy in the 
Yukawa couplings (although such a strong hierarchy already exists in the Standard Model 
Higgs couplings). 

The strong hierarchy in the Yukawa couplings in the allowed solutions implies a very 
specific flavour structure and therefore such models can be tested when more HERA data 
are available. We have considered, as a specific example, the case where the leptoquark is 
dominantly coupled to a positron and a strange squark. 

Note added: Since the HI and ZEUS results were made public, there have been many 
theoretical papers attempting to explain the apparent event excess in terms of leptoquark- 
type particles |^TJ. The present study overlaps with some of these in some respects, but 
to our knowledge our calculation of the HERA cross sections is the first to include the 
0(as) QCD corrections. 
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Appendix 

A.l Kinematics and leading-order cross section 

We consider the process 

e + (p 2 ) + q(pi) - R(q) (Al) 
Pi +P2 = q, s = (pi+p 2 ) 2 = q 2 = m 2 , pj = pj = . (A2) 
The leading order transition amplitude is 

T B = Xv(p 2 ) Pl,r «(pi) , Pl,r = • (A3) 

The spin- and colour-averaged squared amplitude is 



?2 



\ 2 S 



\Tb\ =—, (A4) 



which gives the leading-order cross sections 

Oo{s) = ^-5{s-m 2 ) (A5) 

r 7rA 2 
o"c P = y £0 (s = ys) = -j^r<?(X) > with m2 = xs • ( A6 ) 

A. 2 Virtual corrections 

In the calculation of the real and virtual corrections, we regulate the infra-red and ultra- 
violet divergences using dimensional regularization with d = 4 — 2e. 

A. 2.1 Loop integrals 

We need the loop integrals over two propagators: 

h ^ = I rt fWF^ = Q -^ v ■ ^ ■ < A7 > 

Q_ q2 = i7r */2r(l + e )(-g 2 - , (A8) 



13 



UV,IR 







I 2 (q 2 ;m 2 ) 



euv,m 

•2 (4tt) £ frpUV ,'JJ-i 



+ 2, 



1 A 

y d d k 



if g 2 = , 



1 



(k 2 + irj) [(k - q) 2 + irj\ 2 
1 

{k 2 + [(A; — q) 2 — m 2 + irf\ 



~q»i(q 2 ) , 



him 2 ; m 2 ) 
J 2 (m 2 ; m 2 ) 



g 

^5 m z - 1 



l r 1 



e Jo 

uv 



x — x(l — x) 



m 



,dl 2 (q 2 ;m 2 ^ 



dq 2 



~Qm 2 5 



and the loop integral over three propagators: 



h{q,p) 

h(q,p) 
I 3 (q 2 ;m 2 ) 

him 2 : m 2 ) 



d d k 



{k 2 + irj) \{k — q) 2 + irj] [(k — p) 2 + irf\ 
g 2 ^0,p 2 = 0,(g-p) 2 = 0, 



-Q-? (-\) Ro, Ro = \ ~ % + O(e) , 



q 



d d k 



l 



(k 2 + irj) [{k — q) 2 — m 2 + irj] [(k — p) 2 + irj] 
q 2 ^0,p 2 = 0,(q-p) 2 = 0, 



1 

2^2 



(A9) 
(A10) 
(All) 
(A12) 



(A13) 
(A14) 

(A15) 

(A16) 
(A17) 

(A18) 



It is interesting to compare the massless and massive I3 integrals: the double pole sin- 
gularities have opposite sign and different normalization, and the finite terms are also 
different. We note also the relation 



r(i + e) 



1 



l + e 2 ^)+0(e 3 ). 



r(i-e) 



(A19) 



A. 2. 2 Self-energy corrections 

Self-energy corrections give wave-function renormalization through the relations 



rp 1 T <^S 



2 D dq 2 



1 d£ F 



(A20) 



p 2 =0 
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where Tb denotes the Born amplitude (see eq. (|A3[ )) and — zEs, — «£f are given by the self- 
energy diagrams for massive scalar and massless fermion lines respectively. Introducing 
the auxiliary notation 

-iC F {2K)- d gl^ 



D 1 = (k-q) 2 -m 2 , D 2 = (k- Pl ) 2 , D 3 = k 2 , C 



where Cp = 4/3 is the usual colour factor, we can write 
S s = C I d d k 



D X D 3 



cid d k 2 -^-- D ^ 2 ^ + m2 



Therefore 



and 



where 



-C / d d k 



D 1 D 3 



+ 2C(q 2 + m 2 )I 2 (q 2 ,m 2 



(A21) 



(A22) 



dq 2 
ldS s 



2C 



I 2 {m 2 ; m 2 ) + 2/ 2 '(m 2 ; m 2 ) 



2 dq 2 



CQ r 



CQ, 

as ( 47r/z 2 



rpXJV _ T 

1 o J o 



IK 



r(i + e). 



(A23) 
(A24) 
(A25) 



An \ m 2 

The self-energy of the massive scalar is zero. As in the case of the massless fermion 
correction, the ultraviolet and infra-red divergences exactly cancel. The ultraviolet counter 
term is non-zero 



1 x ct 



a s ( 47r/x 



2\ e 



47r y m 2 

Similarly for the massless fermion line we obtain 
S F = C(2 -d) f d d k 



ra + e) 



- + In — 

e m 



a 



D 2 D 3 
a s (4?r) e /_ T u V 



-2C(d/2-l)-fcI a (pf;0) 



(A26) 



(A27) 



1R 



and therefore 



2 djij 



<*s (47r) e 

'47rrn - 



1 T UV , 1 t^IR 

"2 + 2 



0. 



(A28) 



In dimensional regularization, the self-energy of the massless fermion is zero: the ultravi- 
olet and infra-red divergences exactly cancel. Using the approximate relation ( |A19| ) the 
ultraviolet counter term can be written as 



1 ,CT 

2 hF 



as ( 47r/i 2 



4/T \ m 2 



2e 



(A29) 
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A. 2. 3 Vertex corrections 



The numerator of the vertex correction is given by 

u(p 2 )(-]t + i> 1 )(-]t+2i> 2 )u( Pl ) = u{p 2 )[k 2 -2p 1 k-4p 2 k + 4p 1 p 2 ]u{p 1 ) (A30) 

= T B (k 2 — 4kq + 2p 1 k + 2m 2 ) 
= T B (2D 1 - D 2 + 2m 2 

The vertex correction can then be written as 

2£>i - D 2 + 2m 2 



V 



TuC d a k 



D 1 D 2 D 3 

2m 2 C7 3 (0, m 2 ; m 2 ) + 2CI 2 {p 2 2 ) - CI 2 {m 2 ; m 2 ) 
1 1 



(A31) 



= CQ m ^----^-2 + 0(e) 

Ultraviolet divergences appear in I 2 (p 2 ) and I 2 (m 2 ,m 2 ), and therefore the counter term is 

as ( 4:7i ji 2 ' 



V CT = T B C F 



4tt \ m 2 



r(l + e)(--+ln // 



(A32) 



A. 2. 4 Contribution to the subprocess cross section 

The total contribution of the virtual corrections to the subprocess cross section is 



<7 virt 



= 2a n 



+ + V/T B + + is F + V/T\ 



1„< 1„> 



CT 



(A33) 



= ^C F ^(^)\(l + e)(-\-l-2 + hn^ 



2tx \ m? 
which leads to the final result 

A 2 7r 0:5 / 4nfi 



2e 



virt 



= cv 



4s 2tt 



f^ 2 \ 1 K M 
\ m 2 J T(l - e) q ' 



where 



(A34) 



(A35) 



A. 3 Real contributions 
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A. 3.1 Kinematics, matrix elements, counter terms 

We consider the processes 

e + (p2) + q(pi) -> g(k) + LQ(q), 
e + (p 2 ) + g(pi) ^ q(k) + LQ(q) . 
The kinematical variables are 



(A36) 
(A37) 



q = m , s = 1p\P2 = m + 2kq , t = —2kpi , u = —2kp 2 . (A38) 
The spin- and colour-averaged matrix element squared of process ( [A36Q is given by 

T 2 



A 2 A 2 
\T q \~ = — C F g 2 s Tp(p 2 ,Puk,q) q and a {0) (s) q = — C F g 2 s R [ip(p 2 ,Pu K q) q ) 



(A39) 



where R[ ] denotes the phase-space integral times the flux factor fi 2e /(2s). From explicit 
calculation of the two Feynman diagrams we obtain 



^(P2,Pi;fc,9)« 



m 



m 



2m 2 2kq 



(kq)(kpi) (kq) 2 kp\ 
2m 2 



kpi 



(A40) 



kq 



kq 
kpi 



Similarly, for the process (|A37| ) we obtain 

- A 2 



IT I 
\ ± g\ 

where due to crossing symmetry 



4 



C F g s ip(P2,Pi; k,q) t 



1 - e 



Ip2(p2,-k;-pi,q)c 



(A41) 



(A42) 



Note that here we take into account the fact that the spin average factor for an initial 
gluon is 2(1 — e). The subprocess cross section is obtained by adding the MS counter term 



CT a S 



^(4vr) e -— ^- Cdz P b/a (z,0) af(zs) 
2ir el 1 — e Jo 



(A43) 



where Pb/ a (z,e) denotes the e dependent and spin-independent Altarelli-Parisi splitting 
functions. In the present context we need 

l + z 2 



P q/q( Z , e ) = C F 



P q /g{z,e) 



1 - e 



:i _ z) -<!-») + 2*d -*) 



z 2 + (l- z) 2 - e 



(A44) 
(A45) 
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To perform the phase-space integrals it is convenient to introduce energy-angle variables 
in the subprocess centre-of-mass frame: 



Pi 



0,0,1), p£ = ^(1,0,0,-1) 



(A46) 



k = — (1 — z)(l, 0, sin^, cos 6) , q = pi + p 2 — k . 
2 



Changing variables to 



1 



z = m fs and u — — (1 — cos 



gives, for the d- dimensional phase space integral times flux factor, 

1 /47T/i 2 V 1 rl 



R 



du2u~ e (l-u)- e (l- z 



,l-2e 



32tts \ s 2 J T(l - e) Jo 

In terms of the variables z and u the Lorentz invariant scalar products are 

s s 
m 2 = zs , kpi = — z ) i kq = -(I — z) and k% = s(l — z) 2 u{\ — u) 

A. 3. 2 Contribution to &\ from the e + q — > gLQ process 

In terms of the cms variables, the i/j q function (see eq. (|A40 ) becomes 



(A47) 



(A48) 



(A49) 



ipipi,p 2 ;k,q) q = 4 



l + z 2 
1 - z 



e(l-z) 



1 1 



2u{\ - z 

Introducing notations for the integrals 

J_ a = / duu~ x ~Hl-u) 
Jo 

I = [ duu'Hl -u)- e = 1 + 2e + 
Jo 

and for the functions of z 



(1-z) 2 2 2 



(A50) 



1 7T 2 

e 6 



\o 



:i-z) 



2e V ; (1 - z)+ V 1 ~ z 
Vx = (1 - z) 



then 



1 /4tt^ 



r i 



J_ 1 v_ 1 (l + 2 2 



(A51) 
(A52) 

(A53) 
(A54) 

(A55) 



eI- 1 V 1 -2zV- 1 I -I V 1 



and the unsubtracted cross section becomes 

A 2 7T a,s ( AlX^ 



&lq — Cf 



As 2tt V s 2 J T(l - e) 



K, 



where 



5(1 -z) 
1 I+2 2 



1 1 7T 2 



e 2 e 



e(l-z)+ (I-2) 



22 /ln(l - z) 

+ 2 ' 



1-2 



(1 + 2 2 ). 



The MS counter term cross section (see eqs. ( |A5| . [A43|) ) for this process is 



a 



CT 
1<Z 



A 2 7T a s / An /x 2 



4s 27T 



r(i-e) 



1 1 + z 2 



6(1-2) 



+§5(1 -2) -In ^-^(2, 0) 
Finally for the subtracted partonic cross section we get 

A 2 7T CY S / 4lTfI 2 Y 1 



■sub 



As 2tc \ m 2 J r(l-e) 9 



where 



JT« = 6(1- 



1 5 7T 2 3 , Ut R 

-H h2- — - -ln^ IK 



2e 



6 



2z 



(1-2). 



/ln(l-2)\ (1 + /) _ ln /4 1±4_ 



A. 3. 3 Contribution to di from the e + g — ► gLQ process 

In terms of the cms variables, the ip g function (see eqs. ( |A40| , [X4"2h ) becomes 



ip(p 2 ,Pi] k,q) 



l-e 



where 



2u(l - z) 



2z 2 



ipi(z,u) + ifa(z,u) 



if>i{z,u) = 1-22+- 

1 — u(l — z) 

^(2,0) = 2 2 + (l-2) 2 -e, 



(l-e)u(l-z) 



ip 2 (z,u) 



1 



+ 



2 [i-u(l-z)} 2 l-u(l-z)' 
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The phase-space integral over ip g including the flux factor can be written as 



1 /47T/X 



2\ € 



where 



du 



ins \ s 2 J T(l — e) 



1 



-5(u) + — (1 + e)[l - 2eln(l - z)]Mz,u) 



(A65) 



(A66) 



+2(1 -z) ^ duip 2 {z,u) 



which leads to 



K„ 



z 2 + (1 - zf 



+ 21n(l -z) 



z 2 + (l- z) 2 



(A67) 



+Az - Az 2 + 2z\nz- 2z 2 In 2 z } . 



In this case the counter term simply cancels the 1/e term and introduces the usual In /z 2 

(A68) 
(A69) 



factor of the MS prescription. Finally we get 

Tr 

where 



sub nn X27Ta S 



-K n 



In 



s(l-z) 



As 2tt \ m 2 J T(l -e)" 3 ' 

+ 2*(l-z)(2 + lnz). 



z 2 + (l-z) 2 



A. 3. 4 Finite subprocess cross sections 

The finite subprocess cross sections are obtained by adding the virtual and real con- 
tributions, whereupon the remaining soft and collinear singularities cancel: 

A 2 7r as 



As 2-k 
\ 2 tt as 

U~27T 



C F K q 



T R K g 



(A70) 
(A71) 



where 



K q = K^+kM 



2z An(l-z)> 



l-z) 



1 - z 



1 • ; 2 S-ln / ' l,i ' 



s a-z) 



(A72) 



and Kg is given in eq. ( A69 ). Finally, combining these with the appropriate parton 
distribution functions gives the e + p cross section of eq. fll6|) . 
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A. 4 Calculation of (k^a qj9 

The average leptoquark transverse momentum squared is defined by the weighted subpro 
cess cross sections (see eq. (|A49|) ) 



{k 2 T a q ) 
(k T a g ) 



A 2 



7T 



4 

A 2 7T 



C F g 2 a R[{l-z) 2 u{l-u)^ q {z,u)\, 
T R g 2 s R[{l-z) 2 u{l-u)^ 9 {z,u)}. 



(A73) 



The integrations are finite in four dimensions and can be trivially performed. One obtains 

A 2 7r a s 1 



(k 2 T a q ) 
(k 2 T a q ) 



C F ^- 



4 2tt3 
A 2 7r as 



+ , (A74) 

i " 2 
- - 3z + - 

4 "2vr 

Folding with the parton distribution functions and dividing by the total cross section gives 



1 z 

3z H h 2z 3 - Az 2 In z 

2 2 



the result given in eqs. (|20| , |2T|j22]) . 
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scalar leptoquark production 
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Figure 1: Scalar leptoquark production cross sections in NLO QCD for e + p collisions at 
yfe = 300 GeV. 
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Figure 2: Average scalar leptoquark transverse momentum for e + p collisions at ^/s = 
300 GeV. 
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Figure 3: Scalar and vector leptoquark production at the Tevatron pp collider. 
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